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SPHEROMAK  TILTING  INSTABILITY  IN  CYLINDRICAL  GEOMETRY 

I.  INTRODUCTION  AND  SUMMARY  , 

“  * 

The  compact  torus,  or  taere  speeificallyyspheromak ,  approach  to  mag¬ 
netic  fusion  has  numerous  advantages  as  a  reactor  over  tokamak  designs. 
Most  important  of  these  are  relative  compactness,  the  fact  that  a 
toroidal  blanket  is  not  necessary,  and  the  absence  of  externa]  toroidal 
field  coils. \  The  last  of  these  can  be  expressed  by  saying  that  the 
engineering  beta  (ratio  of  plasma  pressure  to  magnetic  pressure  at  the 

coils)  can  be  quite  high. 

v  1  ?  ' 

Recently ,' theoretical  work  by  Rosenbluth  and  Bussac ’  has  shown  that 

a  spheromak  plasma  in  a  tight  fitting  spheroidal  shell  is  unstable  to  a 

magnetohydrodynamic  mode  which  they  call  the  internal  tilting  mode  if 

the  boundary  is  slightly  elongated,  i.e.  prolate. 

More, recently ,  a  spheromak-type  configuration  has  been  formed  at 

the  University  of  Maryland  by  a  combination  of  Z-pinch  and  theta  pinch 

•  3 

technologies.  Although  the  spherical  wall  equilibrium  of  Ref.  2  des¬ 
cribes  this  device  in  a  qualitative  fashion,  a  more  appropriate  theoret¬ 
ical  model  for  the  Maryland  device  is  an  equilibrium  with  cylindrical 
conducting  walls  and  endplates.  In  addition,  we  feel  that  it  is  of  value 
to  show  that  the  tilting  instability  is  not  specific  to  systems  with 
conducting  walls  that  are  almost  spherical.  Detailed  comparison  of  the 
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model  with  the  experimental  device  and  the  observed  tilting  instability 

I4 

with  the  theory  has  been  presented  elsewhere.  In  this  paper  we  describe 
in  more  detail  two  of  the  theoretical  approaches  used  and  their  justifi¬ 
cation.  v 
/V 

Our  first  approach,  discussed  in  Secs.  II  and  III  is  that  taken  in 
Ref.  2  for  the  spherical  geometry,  namely  the  Wolt  jer-Taylor'’ approach 
of  computing  nonaxisymmetric  force  free  equilibria.  We  show  in  Sec.  II 
that  the  existence  of  nonaxisymmetric  equilibria  with  eigenvalue  less 
than  that  of  the  axisymmetric  equilibrium  is  a  necessary  and  sufficient 
condition  for  instability  if  no  mode  rational  surface  occurs  in  the 
plasma.  Our  geometry,  in  which  the  normal  component  of  the  magnetic 
field  is  required  to  be  zero  on  cylindrical  walls  of  radius  r  =  a  and 
length  z  =  L,  is  inherently  more  difficult  than  that  of  Ref.  2  ,  where 
expansions  can  be  performed  about  an  analytically  tractable  spherical 
boundary.  Nevertheless,  we  are  able  to  compute  nonaxisymmetric  equili¬ 
bria  for  every  value  of  elongation  L/a.  Our  results  show  that  an  in- 
stability  with  m  =  1  (perturbed  quantities  behave  as  e  ,  where  0  is 
the  toroidal  angle)  occurs  for  L/a  >  1.67  ±  1$-  We  identify  this  as  the 
tilting  mode.  Furthermore,  it  appears  that  a  second  m  =  1  instability 
occurs  for  L/a  >  L . lh ,  but  this  mode  appears  to  have  a  smaller  growth 
rate  than  the  tilt.  We  also  find  that  the  m  =  2,  3,  and  *4  equilibria 
as  well  as  those  with  m  *•  <*>  have  higher  energy  than  the  axisymmetric 
equilibrium,  implying  that  all  modes  except  m  =  1  are  stable. 

Our  second  approach,  discussed  in  Sec.  IV,  is  by  means  of  a  linear¬ 
ized  time  dependent  modified  magnetohydrodynamic  code  written 
for  force  free  geometries.  The  results  of  this  code,  specifically  the 
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marginal  stability  point  for  the  m  =  1  tilting  mode,  agree  to  within 
a  percent  with  the  equilibrium  analysis. 


Nonaxisymmetric  force  free  equilibria  for  this  geometry  have  been 

computed  by  a  somewhat  different  method  by  An,  et  al.7  This  paper  and 

8  9 

Ref.  7  were  both  motivated  by  results  of  Marklin  and  of  An.  Marklin 
utilized  a  linearized  time  dependent  magnetohydrodynamic  code  which 
yielded  tilting  instability  for  L/a  >  1.67  and  also  showed  the 
existence  of  an  unstable  m  -  1  internal  kink  mode  when  the  safety 
factor  q  is  greater  than  unity  on  mag  axis.  On  the  other  hand,  the 

9 

work  of  An  showed  that  nonaxisymmetric  force  free  equilibria 
corresponding  to  the  tilting  mode  are  not  obtainable  by  the  method 
of  separation  of  variables.  As  shown  in  Sec.  Ill  and  Ref.  7, 
however,  nonaxisymmetric  solutions  yielding  Marklin' s  tilting  mode 

Q 

result  do,  in  fact  exist,  but  are  not  separable. 


II.  THE  TAYLOP-WOLTJER  APPROACH 

The  Taylor -Wolt jer  approach  is  based  upon  extremizing  the  magnetic 
energy  W  =  J  %  B  dJx  with  K  =  J  A-B  d  x  held  constant.  The  quantity, 

K,  which  is  an  invariant  of  motion  of  ideal  magnetohydrodynamics ,  is 
called  the  magnetic  helicity.  Introducing  a  Lagrange  multiplier  |a  for 
the  constraint,  the  first  variation  of  W  -  |j  K/2  gives 

7  xB  =  uB  (1) 

/V  »v 

where  ju  is  a  constant.  For  this  special  class  of  force  free  equilibria, 
we  have  7xA  =  uA  +  from  which  we  find  W  =  u  K/2.  (We  need  not  con- 

rs*  /V 
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aider  the  complications  that  arise  in  a  system  in  a  toroidal  container 
because  our  system  is  simply  connected.)  Therefore,  the  solution  to  (l) 
with  the  smallest  value  of  n  will  have  the  lowest  energy  and  might  be 
expected  to  be  stable.  We  shall  see  later  that  n  may  be  assumed  to  be 
positive.  We  cannot  expect  this  analysis  to  extend  into  a  vacuum  region, 
where  K  is  not  conserved.  Just  as  obviously,  (l)  cannot  be  satisfied  in 
a  vacuum.  We  therefore  restrict  our  attention  to  internal  modes,  i.e. 
we  assume  that  there  is  a  conducting  wall  at  the  plasma  boundary,  r  =  a 
and  z  =  0,  z  -  L. 

Indeed,  the  second  variation  of  W  -  ^  K/2  is 

6W*  =  J(6B?  -(u6A  •  6B)  d3x.  (2) 

It  is  easily  seen  that  the  usual  potential  energy  6W  of  linearized  ideal 
magnetohydrodynamics11  for  equilibria  satisfying  (l)  and  with  zero  pres¬ 
sure  is  exactly  6W*,  but  with  the  added  condition  that  the  perturbed 
vector  potential  must  be  expressed  in  the  form  of  a  plasma  displacement 
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as  6A  =  xJ5.  Therefore  the  ideal  magnetohydrodynamic  energy  principle 
is  expressed  in  terms  of  6W*  with  the  added  constraint  B*6A  =  0.  This 
implies  6W  a  6W#,  which  guarantees  that  the  condition  6W*  >  0  is  suffi¬ 
cient  for  ideal  magnetohydro dynamic  stability.  Also,  it  is  shown  in  the 
Appendix  that  the  condition  ^B-SA  =  0  can  be  satisfied  by  a  gauge  trans¬ 
formation  except  at  a  mode  rational  surface.  Therefore,  if  no  mode 
rational  surface  occurs  in  the  plasma,  6W*  =  6W  and  we  have  a  necessary 
and  sufficient  condition  for  stability. 

We  shall  see  that  all  solutions  of  (l)  may  be  considered  to  have 

im  0 

only  one  Fourier  component  e1  .  [incidentally,  the  transformation 
m  -•  -m,  Bg  -•  -B0,|u**  *-  H  leaves  (l)  invariant,  proving  that  n  can  be 

assumed  to  be  positive.]  This  property  also  holds  for  the  normal  modes 
in  the  system  whose  fields  are  given  by  the  axisymmetric  (m  =  0)  solu¬ 
tion  to  (l).  Therefore,  if  there  exists  a  solution  6J3  of  (l)  with  eigen¬ 
value  ju,  we  have  6^B  *  £l6A  (plus  a  gradient),  from  which  we  conclude  that 

6W*  =  p.  (Ja  -  n  )  J  6A^  d^x  .  (3) 

Therefore,  if  (J  is  less  than  the  eigenvalue  (i  of  the  axisymmetric  solu¬ 
tion  of  (l),  and  if  no  mode  rational  surface  (i.e.  where  mq  is  an  integer, 
where  q  is  the  usual  safety  factor)  exists  in  the  plasma,  the  axisymmet¬ 
ric  equilibrium  is  unstable.  We  prove  in  the  Appendix  that  the  usual 
boundary  condition  at  a  conducting  wall  at  the  edge  of  the  plasma, 

A 

nxfiA  =  0,  is  easily  satisfied  by  a  gauge  transformation  which  is  compati¬ 
ble  with  the  gauge  in  which  B  •  6  A  =  0. 

Finally  we  can  prove  that  the  solution  of  (l)  with  minimum  eigen¬ 
value  u  is  stable.  If  we  assume  for  contradiction  that  this  minimum 
*o 


energy  state  is  unstable,  it  follows  that  there  exists  some  trial  func- 

2  3 

tion  6A  for  which  6W*  is  negative.  Then,  minimizing  \  -  SWVj'SB  d  x, 

we  find  that  Vx6B  =  ^Q/(l  -  X)6<B,  where  \  <  0  by  assumption.  But 

u  /(l  -  \)  <  u  ,  which  is  a  contradiction, 
o  'o 

Taylor's  celebrated  assertion^  that  the  stability  of  resistive 
modes  can  be  determined  by  computing  the  energy  (i.e.  the  eigenvalue 
pT )  of  nonaxisymmetric  states  is  equivalent  to  the  assertion  that  6W* 
without  constraint  is  the  relevant  energy  principle  for  resistive  modes. 
While  this  is  an  intriguing  proposition,  we  shall  see  that  there  is  no 
need  to  invoke  it,  because  the  tilting  mode  (m  =  l)  becomes  unstable 
for  parameters  with  which  no  mode  rational  surface  exists. 


III.  SOLUTION  BY  SERIES  EXPANSION 

Proceeding  to  obtain  solutions  of  (l),  we  invoke  the  representa- 


B  =  zxtfiji  +  l/jj.  Vxzxtfiji  , 


which  satisfies  (l)  if  \(i  is  a  solution  of  the  scalar  Helmholtz  equation 
2  2 

(V  +  |a  H  =  0.  The  boundary  conditions  =  0  at  r  =  a  and  =  0  at 
z  =  0,  z  *  L,  become 


—  |  =  0 
r  (j  drdz 


at  r  =  a 


at  z  =  0,  z  =  L. 


Proceeding  by  the  method  of  separation  of  variables,  we  find  that  solu¬ 
tions  of  the  differential  equation  with  condition  (5b)  are  of  the  form 

OO 

r  A  T  /,  ,  \  •  imu  ^  a  m  iu  z  t  c\ 

w  =  E  A  J  (k’r)  sm  kz  e  +  A  r  e  e  ^  ,  (6) 

,  n  m  o 

n=I 

2  2  ^ 

where  k'  =  (n  -  k  )  2  and  k  =  nTT/L.  [The  Bessel  function  is  understood 
to  mean  I  (  yk  -  u.  r)  for  k  >  a.]  The  last  term  in  (6)  has  B  identi- 
cally  equal  to  zero,  and  therefore  trivially  satisfies  (5b).  Now, 
imposing  (5a)  we  have  the  condition 


-  —  2  A  J  (k’a)sin  kz  -  —  2  A  kk'J  (k'a)cos  kz 


a  ,  n  m 
n=l 


H  ,  n  m 
n=l 


2imA  am_^  e^z  =0 


for  all  z  between  0  and  L.  We  have  also  started  directly  from  the  three 


>> 


components  of  (l)  and  obtained  (7),  showing  that  the  representation  (b) 
contains  all  solutions  of  (l)  satisfying  the  boundary  conditions  which 
we  impose. 

For  m  =  0  we  find  that  only  one  coefficient  A  is  nonzero  and 

n 


k'a  =  j  ^  ,  the  JL  *th  zero  of  Thus  we  have 

\i  =  (nV/L2  +  J^/a2)15  (8a) 

iji  =  (|a/k')  Jq  (k'r)  sin  (nnz/L),  (8b) 

B  =  -k  J,(k'r)  cos(nn  z/L) ,  (8c  ) 

r  l 

B„  =  n  J.  (k'r)  sin(nnz/L),  (8d) 

U  1 


B  =  k'  J  (k'r)  sin(nrrz/L).  (8e) 

z  o 

The  term  in  (6)  proportional  to  Aq  does  not  contribute  to  V  for  m  =  0. 
The  axi symmetric  equilibrium  with  the  smallest  eigenvalue  n (hence  the 
lowest  energy)  has  JL  =  n  =  1.  The  flux  surfaces  for  this  state  are 
shown  in  Fig.  1.  The  poloidal  flux  is 

rAg  =  rJ1(J11r/a)  sin(nz/L),  (9) 


and  the  safety  factor  q  =  (B-V9 /B-Vr)dr  equals 


q  = 


^rA8 

TTk 


/ 


dr 


2t  ,,  .  v2  2.2 

r  Ji  ^nr'a  -  r  A  q 


1/2  ’ 


(10a) 


8 


t±L 

n2 


/ 


_d£_ 


p2jl(jllP)2 


1/2 


(10b) 


where  and  are  the  two  points  where  the  flux  surface  crosses  the 
midplane  z  =  L/2,  p  =  r/a,  and  =  rAg/a  the  normalized  flux.  From 
(10b)  it  is  easily  seen  that  q  factors  into  a  product  of  a  term  depend¬ 
ing  on  the  normalized  flux  and  another  term  dependent  on  only  the 
elongation  L/a.  That  is,  q(rAg)  is  self  similar  for  every  L/a.  Further¬ 
more,  we  have  shown  that  q  decreases  monotonically  with  distance  from 
the  magnetic  axis  and  that 


qs  =  n/Uk  =  (1  +  2  a2)  A, 


(11a) 


^s^O,! 


qs/0.6. 


(lib) 


where  qg  is  evaluated  at  the  separatrix  rAg  =  0  and  q&  is  evaluated  at 
the  magnetic  axis,  and  jQ  ^  is  the  first  zero  of  Jo-  The  q  profile 
monotonically  increases  with  L/a  (because  the  poloidal  connection  length 
increases)  and  the  q  =  1  surface  first  appears  at  the  magnetic  axis  for 
L/a  =  1.793.  Notice  that  our  model  has  somewhat  more  shear  than  the 
spherical  model  of  Ref.  2. 

For  m/0,  Eq.  (7)  cannot  be  satisfied  by  a  single  n,  because  and 

J'  cannot  both  be  zero  for  the  same  argument.  Similarly,  a  finite  set 
m 

of  nonzero  cannot  satisfy  (7)  because  any  such  finite  set  of  trigono¬ 
metric  functions  are  linearly  independent.  Finally,  we  observe  that  the 
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infinite  set  of  functions  sin(nTTz/L),  cos(nnz/L)  are  not  linearly 

j.  Q  ^ 

independent.  An  obvious  way  to  approach  (7)  is  to  expand  coskz  and  e  M 
in  a  sin  series,  i.e.  extending  all  functions  to  be  odd  on  the  interval 
-L  <  z  <  L.  These  expansions  have  Fourier  coefficients  that  decrease  as 
n~\  We  use  finite  elements  rather  than  sin  functions  as  basis  functions 
to  avoid  these  convergence  problems.  That  is,  we  truncate  the  series 
(7)  at  n  =  N,  producing  a  set  of  equations  with  2N  +  2  real  unknowns 
(since  A^,  A^,  . . . ,  A^  are  complex).  We  then  evaluate  (7)  at  2N  +2 
points  in  the  interval  0  <  z  <  L.  The  resulting  (2N+2)  x  (2N+2)  matrix 
is  a  complicated  function  of  the  eigenvalue  |i .  We  compute  numerically 
the  determinant  d(p )  as  a  function  ofju.  This  determinant  exhibits 
nonphysical  roots  d(  n)  ~  |  n  -nrr  / L |  .  This  is  because,  as  |U~*n n  /L, 

the  e1^  Z  solution  of  (6)  is  redundant  [leading  to  two  pairs  of  identi- 

m  /O 

cal  columns  in  the  matrix]  and  because  J  (k'a)  ~  k '  J '  (k  'a)~(|j-nTT  /L) 

m  m 

[there  are  two  columns  proportional  to  this  factor.]  The  other  roots  of 
d(n )  are  double  roots,  due  to  the  symmetry  B-*©  +  n/2m  [which  exchanges 
real  and  imaginary  parts  of  A^,  A-^,...,Ajj.]  The  determinant  d(n  )  does 

not  reach  zero  for  finite  N.  For  N  =  10,d(|j  )  exhibits  fairly  sharp 

-2 

dips  down  to  10  compared  to  neighboring  values.  For  N  =  15 ,  the  dips 

-b 

are  more  pronounced,  around  10  .  This  shows  that  quite  a  few  Fourier 

harmonics  are  necessary  to  describe  the  solution  accurately.  The  com¬ 
plexity  of  these  eigenfunctions,  in  particular  the  fact  that  they  are 
not  separable,  is  related  to  the  fact  that  (5)  represents  an  unconven¬ 
tional  set  of  boundary  conditions  for  an  elliptic  equation. 

In  Fig.  2  we  show  the  first  two  physical  roots  for  m  =  1  and  the 
first  root  (8a)  for  m  =  0,  as  a  function  of  the  elongation  L/a.  There 
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are  crossings  at  L/a  =  1.67  and  L/a  *  **.lk.  Since  no  mode  rational  sur¬ 
face  q  =  n  exists  for  L/a  <  1.793  (the  maximum  q  at  L/a  =  1.67  is  0.944) 
this  implies  that  all  m  =  1  modes  are  stable  for  0  <  L/a  <  1.6?  and  that 
an  m  =  1  instability  exists  for  1.67  <  L/a  <  1.793.  The  second  crossing, 
which  occurs  when  the  q  =  2  surface  is  present  (q  has  the  range  l.?9  < 
q  <  2.1 h)  indicates  but  does  not  prove  the  existence  of  another  m  =  1 
instability.  There  does  not  seem  to  be  a  further  crossing  of  the  two 
m  =  1  states,  and  therefore  we  believe  that  the  first  m  =  1  mode  is 

always  dominant  in  terms  of  growth  rate  and  saturation  amplitude.  Be- 

2 

cause  this  mode  has  similar  structure  to  that  of  Rosenbluth  and  Bussac 
in  spherical  geometry,  and  because  it  becomes  unstable  at  a  similar 
elongation  and  when  no  mode  rational  surface  exists,  we  refer  to  this 
as  the  tilting  mode. 

The  lowest  eigenvalues  for  the  m  =  2,  3,  and  l*  modes,  as  well  as 

for  m  =  0  are  shown  in  Fig.  3.  Since  there  is  no  crossing  with  the 

lowest  m  =  0  eigenvalue,  it  follows  that  there  are  no  m  =  2,  3  or  i 
instabilities.  [Note  that,  as  L/a-*  0,  the  second  term  of  (7)  becomes 
dominant,  indicating  that  the  eigenvalue^  should  become  independent 
of  m  in  this  limit.  The  results  of  Fig.  2  and  3  show  that  this  is  in¬ 
deed  the  case.]  Furthermore,  for  m  00  and  Aq  =  0,  (7)  is  satisfied 

with  a  single  nonzero  A  ,  giving  |U  =  (n  n  /L  +  J^/a  )  •  The  lowest 
energy  state  has  n  =  X  =  1  .  The  m  =  U  results  of  Fig.  3  agree  with 
this  value  within  a  few  percent  and  we  therefore  conclude  that  the 
axisymmetric  equilibrium  (8)  is  unstable  only  to  m  =  1  modes. 


IV.  LINEARIZED  CODE 


Here  we  present  an  alternate  approach  to  that  of  Sec.  III.  First, 
notice  that  the  linearized  equation  of  motion  of  ideal  magnet ohydrody- 
namic  for  force  free  equilibria 

p  d^/dt2  =  6^j  x  B  +  ^  x  6B  (12) 

takes  the  form 

(p/B2)  (d2/dt2)6A  =  (13) 


for  equilibria  satisfying  (l),  where  6A  =  §  x  B  and  ±  signifies  the 
components  perpendicular  to  the  equilibrium  field  B.  Thus  we  may  con- 
sider  6A  and  not  §  our  basic  dependent  variable;  clearly  B-6A  =  0  for 


all  time  if  B*6A  and  (d/dt)B-6A  are  both  zero  at  t  =  0.  This  set  of 

A#  /V  /V 

equations,  of  course,  has  a  variational  principle  with  6W*,  with 
normalization  T  =  J  d^  x  p|6A|2/B2  and  with  constraint  B-6A  =  0.  The 
modified  set  of  equations 


(p/B2)(d2/dt2)  6 A  =  pi6B  -  6^ 


(HO 


has  the  same  variational  principle  but  without  the  constraint.  As  dis¬ 
cussed  in  Sec.  II,  the  stability  criterion  obtained  by  integrating  (13) 
or  ( lU )  must  be  identical  if  no  mode  rational  surface  exists  in  the 
plasma.  Furthermore  ( lb )  has  the  advantage  of  producing  nonaxi symmetric 
solutions  to  (lU)  at  marginal  stability.  [Taylor's  assertion^  mentioned 
earlier,  is  therefore  equivalent  to  the  statement  that  tearing  modes  are 


12 


stable  in  an  axisymmetric  system  satisfying  (1)  if  and  only  if  (14)  is 

stable,  regardless  of  the  existence  of  mode  rational  surfaces.] 

Our  code  integrates  (14)  in  time  until  one  unstable  mode  if  any 

dominates  all  others.  Runs  are  terminated  when  the  growth  rate,  based 

upon  the  kinetic  energy,  varies  by  less  than  a  percent.  Results  from 

2 

a  series  of  such  runs  with  normalization  p/B  =  const,  are  shown  in 

Fig.  4,  for  m  =  1.  The  marginal  stability  point  at  L/a  1.65  agrees 

well  with  that  shown  on  Fig.  2.  The  perturbed  magnetic  field  of  this 

unstable  mode  near  marginal  stability  (L/a  =  1.75)  is  shown  in  Fig.  5. 

Notice  that  the  real  part  of  (i.e.  the  cosm©  component)  looks 

like  a  pure  sine  wave  as  a  function  of  z,  but  with  a  period  slightly 

less  than  L,  whereas  the  imaginary  part  of  6 (i.e.  the  -sinm 

component)  seems  to  be  a  combination  of  sin(irz/L)  and  cos  (2nz/L)  . 

The  other  components  of  also  show  a  complex  structure  of  Fourier 

harmonics  in  z,  including  aperiodic  behavior.  See  especially  the 

imaginary  part  of  SB. ,  which  seems  to  have  a  period  1.16  L.  This  is 

y 

to  be  expected  since  the  m  -  1  solutions  to  (7)  contain  many  harmonics 

iy  z 

as  well  as  the  anharmonic  component  proportional  to  e 
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APPENDIX 


Our  purpose  here  is  to  prove  the  two  propositions  of  Sec.  II.  The 
first  is  that  a  gauge  transformation  can  be  found  so  that  6 A  is  normal, 

A 

i.e.  n  x  5 A  =  0  on  the  boundary,  if  the  volume  is  simply  connected  and 

A 

if  n-6ji  =  0  on  the  boundary.  The  second  proposition  is  that  the  con¬ 
straint  B*6A  =  0  can  be  satisfied  throughout  the  plasma  if  no  mode 

ew  »v 

rational  surface  exists  in  the  plasma  (as  is  the  case  in  the  neighbor¬ 
hood  of  marginal  stability  for  the  tilting  mode). 

To  prove  the  first  proposition,  we  ask  whether  there  is  a  scalar 

A  A 

function  X  such  that  t-Vx  *  -t*6A  =  0  on  the  boundary.  Assuming  X=  0 
at  some  reference  point  on  the  boundary,  we  find 


x<*> 


-(% 

J  x 

~o 


61  . 


The  function  \  so  defined  is  single  valued  because  the  integral  over  a 
closed  path  equals  the  magnetic  flux  through  that  part  of  the  boundary 
encircled  by  the  path. 

To  prove  the  second  proposition,  we  must  find  another  scalar  func- 

A 

tion  X  that  satisfies  B'^x  =  -B-6A  (where  n  x  6A  =  0  is  now  assumed  at 
the  boundary).  Writing  =  Vijf  x  V8  +  g(ij/  )V0  is  now  the  poloidal  fiux 
rA^iwe  find 

J  t ^  It  - -2'^  •  (#-1) 


where  cp  is  a  coordinate  along  the  flux  surfaces  i|f  =  const.,  and  J  is  the 

-1  2 

Jacobian  (V^  xVQ  -Vcp)  .  We  can  always  assume  J  =  r  K(| ) ,  where 

p 

K( + )  *  (l/2  n)  $  d£/(r  B  )  and  B^  is  the  magnitude  of  the  poloidal  field. 
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,  we  can  assume 


2  _  Q 

and  that  cp  varies  from  0  to  2n  .  Also,  since  6A~  e 

/v 

im  0 

X  *  is  also  possible  to  add  an  arbitrary  function  of 

$.)  Substituting  in  (A.l),  we  find 

dx/dcp  +  imq(>f)x  =  -r2K(it  )B-6A,  (A. 2) 

where  q(i|i )  =  9(^)K(i})  is  the  safety  factor  (l/2rr  )  §  (B-V  0  )dcp /(B^Vcp  ). 
Expanding  X  and  the  right  hand  side  of  (A. 2)  as 

2  Xn(*)e'ln;p 
n=-°° 

and 

Z  f  (Y)e~incp 
n—  ° 

and  substituting  in  (A. 2),  ve  find 

Xn(^  =  1  fn(t)/(n  ~  Diq(iji)).  (A. 3) 

If  ♦  *  ♦q  is  a  mode  rational  surface,  i.e.  has  n  -  mq(i|io)  =  0  for  some 

n,  we  must  have  f n ( *j> Q )  =  0.  If  no  mode  rational  surface  exists,  on  the 

other  hand,  there  is  no  constraint  on  6A,  i.e.  on  the  coefficients  f  (i|i). 

**  n 

It  is  obvious  that  this  second  gauge  transformation  does  not  destroy  the 
property  guaranteed  by  the  first.  We  wish  to  emphasize  that  these  two 
properties  are  general  and  do  not  depend  upon  B  being  a  solution  of  (l). 
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L/a 

Fig.  2  —  Normalized  eigenvalue  j/L/ff  as  a  function  of  elongation  L/a  for  the  lowest 
m  ”  0  equilibrium  and  for  the  two  lowest  m  ■  1  equilibria.  Note  the  crossings  at 
L/a  “  1.67  and  at  L/a  ■  4.14 
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L/a 

Fig.  3  —  Normalized  eigenvalue  fiL/ir  as  a  function  of  L/a  for  m  ■  0,2,3,  and  4  equilibria 
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Fig.  4  —  Normalized  growth  rate  y  =  IojIL/V^  for  the  most  unstable  mode  of 
the  system  (14)  with  m  =  1,  as  a  function  of  elongation  L/a,  showing  marginal 
stability  at  L/a  =*  1.65 
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Fig.  5  —  Magnetic  fields  produced  by  the  m  =  1  tilting  instability  from  (14)  near 
marginal  stability  (L/a  =  1.75).  Respectively,  these  are  the  real  and  imaginary  parts 
of  6Br  (a,b),  of  5Bfl  (c,d)  and  of  5BZ  (e,f).  For  the  solid  curve  the  horizontal  axis 
represents  r  (at  z  =  L/2)  and  the  zero  point  is  marked  on  the  left  hand  vertical  axis. 
For  the  dashed  curve  the  horizontal  axis  represents  z  (at  r  =  a/2)  and  the  zero  point 
is  marked  on  the  right  hand  vertical  axis.  The  scale  of  SB  is  arbitrary. 
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